Composite-Fermion Theory for Pseudogap, Fermi Arc, Hole Pocket, and 
Non-Fermi-Liquid of Underdoped Cuprate Superconductors 
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We propose that an extension of the exciton concept to doped Mott insulators offers a fruitful 
insight into challenging issues of the copper oxide superconductors. In our extension, new fermionic 
excitations called cofermions emerge in conjunction to generalized excitons. The cofermions hy- 
bridize with conventional quasiparticles. Then a hybridization gap opens, and is identified as the 
pseudogap observed in the underdoped cuprates. The resultant Fermi-surface reconstruction nat- 
urally explains a number of unusual properties of the underdoped cuprates, such as the Fermi arc 
and/or pocket formation. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.30.+h, 74.72.-h 



Since the discovery of cuprate superconductors, the na- 
ture of low-energy electronic excitations evolving in their 
normal metallic phase has attracted much attention as 
one of the central issues in condensed matter physics. 
One reason for the interest lies in its connection to the 
origin of the high temperature superconductivity itself. 

Electronic states in the underdoped cuprates are un- 
conventional. For example, spin and charge excitations 
are unexpectedly suppressed as "pseudogap phenomena" 
in the normal state. Recent improvement of experi- 
mental tools, such as angle-resolved photoemission spec- 
troscopy (ARPES), has further enabled resolving strong 
momentum dependence of quasiparticles 1 - - — . In partic- 
ular, quasiparticles are hardly observed around antin- 
odal points (tt,0) and (0,7r) in the 2D Brillouin zone 
for the CuC>2 plane. It looks a truncation of a large 
Fermi surface observed in the overdoped cuprates, and is 
sometimes called the "Fermi arc" . More fundamentally, 
the normal state of the cuprates remains a challenge as 
Mott physics in the proximity to the Mott insulator—. 
Although the doped Mott insulators in two dimensions 
have been studied for a long time from various theoret- 
ical approaches^—, the nature of the electronic states 
is not yet fully understood. Recently revealed pseudo- 
gap and arc or pocket of the Fermi surface^— require a 
conceptually deeper understandig for Mott physics. 

In this letter, we elucidate a key role of exciton-likc 
physics on this issue. Excitons are known to be a key 
concept in physics of semiconductor o 1 7 i 18 . The excitonic 
state also emerges in the Mott insulator, for instance as 
the charge transfer excitation at the optical gap edge 
in the cuprates 1 ^— , due to a strong binding of empty 
(holon) and doubly occupied (doublon) sites in half-filled 
Mott insulators. In the doped Mott insulators, in spite of 
screening by doped carriers, the remnant of the binding 
may still remain as weak binding between a doped holon 
and the preexisting doublons similarly to excitons. When 
an electron or hole is added to the doped Mott insula- 
tors, it may appear as a normal quasiparticle extended in 
space. However, an electron (a hole) can alternatively be 
added locally to a holon (doublon) site with a small cost 
of the on-site Coulomb repulsion. This generates a bound 



composite particle (cofermion) consisting of the preexist- 
ing holon (doublon) and the added electron (hole). We 
call such cofermions holo- electron (doublo-hole). 

A cofermion (a holo-electron or a doublo-hole) dynam- 
ically breaks up into (and recombines from) a conven- 
tional quasiparticle and a charge boson. This dynam- 
ical process is naturally interpreted as a hybridization 
between the cofermion and the quasiparticle. Here, we 
show that the resultant hybridization gap offers a nat- 
ural understanding of a number of key properties of 
the underdoped cuprates 4 - such as Fermi pocket or arc 
formation^— , psudogap behavior seen in the single parti- 
cle spectra^ 2 -, specific hea t 23 : 24 , the asymmetric density 
of states (DOS) 2 ^, and violation of the Wiedemann- Franz 
la w 26 ' 27 , without any symmetry breaking. We specifically 
predict that the pseudogap opens as a s- wave-like gap in 
the unoccupied part above the Fermi level contrary to the 
widely assumed d-wave structure. 

We study the Hubbard hamiltonian on a square lattice, 



(1) 



where c i<7 (c ia ) is spin-a creation (annihilation) operator 



at a site i, while h i 



-cLc 4 . 



For the hopping ty, we 



take only —t for the nearest-neighbor and t' for the next- 
nearest-neighbor pairs. 

We first employ the Kotliar-Ruckenstein slave-boson 
formalism 2 ^, while the local Hilbert space of the Hub- 
bard model is expanded not by the original electron c a 
but instead by introducing a fermion f a , which stands 
for the cr-spin quasiparticle, following Refl29l and one 
slave boson for each Fock state as e for the empty state 
(holon) |0), p a for the singly occupied state \cr) (<r=t, 
or I), and d for the doubly occupied state (doublon) 
|tl)- After the mapping, the Coulomb repulsion U is 
now interpreted as a "chemical potential" for d\, while 
the correlation now appears, as we describe below, in 
hopping process of j ia disturbed by slave-boson mo- 
tion under the local constraints ete,-+^ p\ J 'p ia -\-d} i d i =Y 
and fi a ha = p\aPi<j Jr ^ L \^i imposed to keep consistency be- 
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tween the boson and fcrmion Hilbcrt space. In the en- 
larged Hilbert space, these two constraints are assured 

respectively by the Lagrange multipliers A^ and A- 2 ? in 
the Lagrangian as 

£ = E flMfiitii + CijAr^Ar) + 4°\ (2) 
where D^d T - fi+X\J, Qj< T { r )=kA T )ZjA T )i and 

+ fe)^ ( -4iw - ^e 4 )- 1/2 

and s{ 2) = (l - pLAct - dj^i) _1/2 ' following the 
literature^). A part of the Lagrangian contains 
Xf' , X^J and quadratic terms of bosonic fields only2£ as, 
4 0) =E, {e\{r)[d r + Af } ]e, (r) + £ CT pj»[d r + A« - 
+ d\(r)[d T + U + A« - ZrX^M (r)}. To 
take into account the Gaussian fluctuations of the bosonic 
fields beyond the mean-field level^, the Bogoliubov pre- 
scription is useful, where the boson operators 6 4 (b = e,d 
or p CT ) are divided into condensate 60 and fluctuating 
components b i as b\=b +b\, ^=60+6-. 

In this letter, we further make a progress by consid- 
ering low-energy dynamics arising from coupled bosons 
and fermions. First, we reexamine the strong couling 
(U/t — > +00) limit&i£, where adding a er-spin electron is 
possible only to |0), namely, a holon site, to avoid creat- 
ing |ti) (doublon) with the cost of U. Creation operators 
for the electron at |0) are given by composite fcrmion op- 
erators e i f- a . 

When t/U becomes nonzero, an added electron may 
become a coherent quasiparticle. However, we still have 
rather localized character of holons, and it allows alter- 
natively forming a collective excitation of a hole and the 
added electron similarly to an exciton. This collective 
character is clearly distinguished from the conventional 
quasiparticle. In fact this composite fcrmion does not 
have charge in contrast to the quasiparticle. Our crucial 
step is to include dynamics of this composite fermion ex- 
pressed by eJl (j^ej. 

When we impose the local constraints more strictly 
for fluctuating bosons beyond the mean-field level, it 
turns out that, in the hopping process of f a expressed 
by fiJijakjfja, the coefficient £fg =q\ a 9l a + 
3\Pia)(^jPja ^Pja^j) I s dominating (see Appendix. [Taj. 
Here we employ gia={l - P% - el)~ 1/2 and - 

Poct — °W following Refill This vertex stands for 
the backflow consisting of bosons, originating from the 
quasiparticle motions. 

Then we treat coupling of charge bosons and quasi- 
particles in fJ a Cijlfja- which represents a part of the 
electron-electron interactions in the hamiltonian (JTJ), by 
interpreting the form such as fj^b^jfj^ as the decou- 
pled product of C ia =(e i: d^)fl and C %a ^f la {e\,d i ) T . 
Namely, this boson-quasiparticle interaction is equiva- 
lently treated by introducing integrals over the grass- 




FIG. 1: Diagrams for the Dyson equations. Solid and dashed 
lines with arrows represent propagators of the quasiparti- 
cles and cofermions, respectively. Wavy lines stand for the 
charge bosons, and solid lines are the spin bosons. Condensa- 
tions of bosons are represented by lines terminated at crosses. 
Coupling constant g\ a g\otij ls represented by open polygons. 
Here, we do not distinguish holons and doublons. Spins and 
"flavors" of cofermions (ip, x) are also not distinguished in the 
diagram, for simplicity. Filled circles stand for the amplitude 
of the hybridization between quasiparticles and cofermions. 

manian Stratonovich-Hubbard fields T^ .=(V'] ., %L) anc l 
"^'i(r = (V'i(7! Xia) T where the interaction is formally trans- 
formed to the hybridization of T and C as C , n C r . — > 
C ia T jrT + T ia C ja - T ia .Yj CT (sec Appendix. ITE) . 

The newly introduced Grassmann fields ipj a and 
X ia arc physically interpreted as cofermions, the 
holo-electron and doublo-hole, respectively. Then, 

C^Tjct+T^Cjo. is naturally interpreted as breakup 
and recombination processes of the cofermions. After in- 
tegrating fluctuating bosons out, as is mentioned below, 
it results in the hybridization between cofermions \ a 
and quasiparticles / CT . 

We treat Gaussian fluctuations of bosons, and the dy- 
namical coupling between quasiparticles and cofermions 
by using a set of the Dyson equations up to the second 
order of tij, as is depicted in FigfTJ Thick lines (thick 
wavy lines) stand for the dressed Green's functions of the 
charge bosons A ab (r, t)=- (T (t) (0)) (spin bosons 
C oi, (r,r)=-(T^(r)^ t (0))), where a, 6=1, 2, r=i-j, 

(A 1 ,/3 l 2 )=(e 4 ,d l ), and , P%). Here we ne- 

glect the coupling between charge and spin bosons such 
as (pj^e^), which vanishes in Mott insulators and gives 
higher order terms scaled by the doping rate x for |x|<Cl. 
Thick lines with arrows represent the dressed quasiparti- 
cles Q a {r, r). Thin lines (thin wavy lines) represent bare 
propagators of the charge bosons AQ b (r, r) (spin bosons 
Cg fc (r, t)), determined by £g' . Thin lines with arrows 
are bare propagators of the quasiparticles Q Q(T (r, r) de- 
termined by C =J2i 3 /L( r )[A<% + Coatij}f ja (T), where 

Coa=gf a gh(Poaeo + ~d~oPoa) 2 - Tne Lagrangian C is ob- 
tained by decoupling the fluctuating bosons from the La- 
grangian L - (see Eq.©). Thick and thin dashed 
lines with arrows are the cofermion propagators J- cd 
(c,d=ip,x) and F^^Sd/t {(- — » 0), respectively. This 
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FIG. 2: (color online). Calculated spectral functions for 
t'/i=0.25, U/t=12 and £=0.05. (a) Spectral function ,40,0;) 
along lines running from (0, 0), (tv, it) and (n, 0) to (0, 0). The 
dashed line is the Fermi level. We use a finite broadening fac- 
tor 5=0. 05t. (b) Band dispersion of quasiparticle for £=0.05. 
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FIG. 3: (color online). Reconstructed Fermi surfaces and 
spectral functions Af(k,ui) for i'/t=0.25 and U/t=12. (a)- 
(c) Af(k,u))(= — Im[G^/' (k, u})]/tv) at uj = i8. Here we take 
the broadening factor S = 0.05t. Solid and dashed lines il- 
lustrate the poles of quasiparticles. (d) Doping dependence 
of Fermi-surface topology in our theory, (e) Af(k,ui) along 
the symmetry line running from (0, 0) through (n, 0) to (n, n) 
for x — 0.10. Thin solid and dashed white lines illustrate 
poles of quasiparticles. We define Apg as the gap between 
(X (w = 0) and the maximum of the quasiparticle dispersion 
below n along this symmetry line. 



peculiar divergence of Jq is because the Lagrangian docs 
not include cofermions if we neglect the interactions be- 
tween bosons and fcrmions. By solving the Dyson equa- 
tions, we obtain the propagators for the quasiparticles 
and cofermions. Here the bosonic degrees of freedom are 
taken into account in a self-consistent fashion, through 
the cofermion self-energy and the amplitude of hybridiza- 
tion between the quasiparticles and cofermions. 

Now we show how our self-consistent solution of cou- 
pled quasiparticles, bosons and cofermions predicts nor- 
mal state properties. We show the result at U=12t and 
t'=0.25t to get insight into the cuprate superconductors 
by restricting to paramagnetic solutions at temperature 
T=0. First, we give the spectral functions calculated 
from the electron Green's function G a (k,ui) (see Ap- 
pendix. 1 1 c[) ; A(k, w)=-Im [G a (k,tu)] /n. In Fig|3h, we 
show A(k,u) for the hole concentration x = 0.05. Two 
main features are found; the coherent band arising from 
the quasiparticle around the Fermi level and the remnant 
of the upper (lower) Hubbard band at tu>6t (lj<— 2t) 

generated by dynamics of e and rf 30 i 31 . 

Here, we focus on reconstructions of the Fermi surface 
in the coherent band. The quasiparticle Green's function 
is given as 

GW{k,u)~ {lu - Co*e k + n -Zfik^)}" 1 , (3) 

where £/=A(fc) 2 /(7fcu; — a k ) is the quasiparticle self- 
energy arising from the quasiparticle-cofermion hy- 
bridization A(fc). Here the cofermion propagator (7^ — 
ojfc)" 1 is obtained from the expansion around the 
cofermion pole (sec Appendix. 1 1 c[) . The quasiparticles 
Green's function has a hybridization gap due to the hy- 



bridization with the cofermions and Gj shows the diver- 
gence of the quasiparticle self-energy given by S/ at the 
zero surfac o 11 ' 13 " — defined by 7fcu; — afc=0. In our theory, 
the zero surface splits the band dispersion and generates 
a distinct s-wave-like gap (as is seen in Figj^b and sup- 
ported by recent numerical observatio n 1 ^ 15 ) . For small 
hole doping such as a;=0.05, our theory predicts that the 
reconstructed Fermi surface becomes a small pocket as 
we see in Fig(3Ji-[5t. 

Actually, the topological transitions occur at x~0.13 
and x~0.18 (see FigJ5U). Below x~0.13, the Fermi sur- 
face consists of small hole pockets. It is difficult to dis- 
tinguish the pockets from the arc structure as we see in 
Fig|3k. This is because the zero surface near the outer 
part partially destroys the quasiparticles. For 0.13<x, 
large Fermi surfaces appear, instead of Fermi pockets. 
For 0.13<a;<0.18, a hole-like surface centered at (7r,7r), 
and an electron-like one centered at (tt, tt) coexist (see 
Fig©})- However, the electron- like surface is hardly seen 
again because of the nearby zero surface. Our result 
shows a nontrivial topological transition at x~0.13 as 
a consequence of the hybridization with cofermions. 

A gap A PG measured from the Fermi level \x emerges 
near the antinodc, corresponding to the pseudogap in 
the ARPES as we identify in FigEJ 3 - The pseudogap 
Apg is determined by the hybridization gap A(fc), basi- 
cally scaled by a fraction of t, consistently with numeri- 
cal studies^. The doping dependence of Apg is given in 
Figgk, in agreement with the ARPES for LSCC&2 2 .. 

For x<0.13, the density of states (DOS) of the elec- 
trons c ka at the Fermi level, p-p— — Co / d 2 kAf(k, 0)/47r 3 , 
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FIG. 4: (color online), (a) Pseudogap Apg as a function of 
hole doping rate x. Thick solid line is Apg calculated by our 
theory. Closed (red) squares show Apg estimated from the 
ARPES^ for La 2 _ I Sr a; Cu04 (LSCO). (b) DOS of electrons 
vs. x in the present theory (bold solid curve). The thin 
dashed curve stands for the DOS of the non-interacting case. 
Closed (blue) circles are obtained by a linear extrapolation of 
low-temperature normal-state 7 of LSCO in Ref[23| to T— 5-0, 
representing the presumable lower limit. Closed (red) squares 
stand for the estimate for LSCO obtained in Refl24l by the Zn 
doping, which may be an upper limit. All of the present 
results are obtained at t=0.4eV, i'/t=0.25 and U/t— 12. 



optical conductivity^ indeed supports our prediction. 

Our finding is that hidden cofermionic particles called 
holo- electrons and doublo-holes play a key role: The 
cofermions hybridize with the quasiparticles and cause 
a hybridization gap identified as the pseudogap. A num- 
ber of resultant properties consistent with the unusual 
normal states of the cuprates support relevance of our 
cofcrmion theory to physics of the cuprates. 

The authors thank Yukitoshi Motome and Shiro Sakai 
for useful discussions. Y.Y. is supported by the Japan 
Society for the Promotion of Science. 



Appendix 

In this Appendix, we show theoretical details on in- 
troducing the cofermions, construction of quasiparticle 
Green's functions, and supplementary results on the 
quasiparticle density of states. 



1. Details of theory 



is clearly suppressed, as is illustrated in FigHJ}. We 
compare pp with the specific heat coefficient 7 mea- 
sured for LSC O 23 ' 24 by using the conventional relation 
7=7t 2 /9f/3 at T=0. Our 7 is consistent with experi- 
ments. The w-dependence of the DOS shows significant 
asymmetry around w=0 in contrast to the DOS for the 
non- interacting case (see Appendix. [2j ■ This asymmetry 
of the DOS naturally explains the asymmetric tunneling 
spectra observed in the hole-undcrdoped cuptratea 2 ^ (see 
Refsl32l and l33l for different interpretations). 

The present result slightly depends on the choice of 
the parameters. For instance, Apg decreases from the 
present result by an amount ~0.05t at i'/t=0.25 and 
U/t=15 or i'/£=0.15 and U/t=12, while the qualitative 
features are robust. 

Here we note the difference of the cofcrmion from the 
spinoiA Although the cofermion carries a spin but no 
charge as in the spinon, cofermions coexist with quasi- 
particles in different part of energy-momentum space as 
the electron differentiation in contrast to the spinons. 

The present cofcrmion contributes to the entropy and 
the thermal conductivity n in addition to the quasipar- 
ticle. On the other hand, the electric conductivity a is 
contributed only from the quasiparticle. Therefore we ex- 
pect a serious breakdown of the Wiedeman- Franz la w 26 ' 27 
that predicts a universal constant Lo=ir 2 k 2 i /3e 2 for the 
ratio L=k/Tct. Our theory predicts L > Lq. 

We propose to test our specific prediction of the s- 
wave-like gap structure in unoccupied spectra, for ex- 
ample, by improving the low-energy electron spectro- 
scopies, such as the inverse photoemission, the low- 
energy electron diffraction spectroscopy, resonant inelas- 
tic X-ray spectroscopy, or time resolved photoemission 
spectroscopy. The mid-infrared peak and long tail of the 



As our theoretical starting point, we employ the 
Kotliar-Ruckcnstcin slave-boson formalism^ 8 - for the 
Hubbard model, where the local Hilbert space of the 
Hubbard model is expanded by introducing one slave bo- 
son for each Fock state as e for the empty state (holon) 
|0), pa for the singly occupied state \a) (<r=t,or I), and d 
for the doubly occupied state (doublon) |f4). In addition 
to these bosons b (b = e,p„ or d), fermion operator f a 
is introduced to stand for the a-spin QP. The mapping 
between the original electrons and / combined with b is 
given by c\ a = z irT fJ a , where z ia is define d 28 ' 31 as 

with 

5^ = (l-pU^-eK)- 1/2 , (A.2) 

and 

9$ = ^-vlPi<,-$idi)- X/ \ (A.3) 

following the literature^ 8 .. 

We need to impose local constraints to eliminate un- 
physical states. First, only one boson should occupy each 
local state as 

+ Y, pIp™ + = 1. (A.4) 

a 

Second, the number operator of fa- is necessarily given as 

flfi*=plPi* + <%di- (A.5) 

These constraints are represented by integrals over La- 
grange multipliers A| , and , in the path integral of 
the Lagrangian discussed below. 



In the expanded Hilbcrt space, the Lagrangian of the 
Hubbard Hamiltonian is mapped to 

L = E fU T )& 5 v + 4>)M/,>) + 4°\ (A.6) 



where A = d T - fi + \\ 2 J and QjA T ) = hA T )^ 3 A T )- 

contains \\ , and quadratic terms of bosonic 
fields only2& as, 

+ dJ(r)[6V + t/ + Af > - ^ ^]^(r)}.(A.7) 

The on-site Coulomb [/ is now interpreted as a "chemical 
potential" for d\, while the correlation now appears in 
hopping process of j\ a disturbed by slave boson motion. 

a. Bosonic fluctuations 

To take into account the Gaussian fluctuations of the 
bosonic fields beyond the mean-field levels, the Bogoli- 
ubov prescription is useful, where the boson operators 
are divided into condensate components bo and fluctu- 
ating components b i as b\ = bo + b\, b i = bo + b { with 
b = e,d or p a . 

When we impose the local constraints more strictly for 
fluctuating bosons beyond the mean-field level, it turns 
out that the term 



(b-1) - „ 



(1) 

ija 



<?L <?L (pL ^ + ^Piv ) (e]Pj<T + P]a d j ) ( A 



represented by the diagram in Fig|3^ is dominating 
among all the possible diagrams for Qj a in Eq. (IA.6l) . 
Here we employ g\ a = (1 - p 2 - - e^)' 1 and g\ a = 

(1 -via - dlr 1 by following ReflH. 

To elucidate why we retain CiiV' we classify the dia- 
grams illustrated in Fig JS] into four types, categorized by 
time dependence (or frequency dependence) of quasipar- 
ticles and the local conservation of the boson densities: 

(T-l) Diagrams containing external propagators of 
quasiparticlcs, in addition to bosonic propagators 
violating the local conservation of boson densities 
(Figures [SbfSklSH, and [5h). Here the violation 
means that before and after the interactions (rep- 
resented by hexagons), the number of bosons ex- 
pressed by external boson propagators is not the 
same. 

(T-2) Diagrams that contain time dependence of quasi- 
particles, but that do not violate the local conser- 
vation (Figures [5^ andJSJ;). 



:-l) 

1m 



(d-2) 




(h) 



FIG. 5: (a)-(e) Diagrams representing terms in fia-Ojlja- (f) _ 
(h), Examples of the diagrams for perturbation generated 
from the same term as (a)-(e). Solid lines with arrows repre- 
sent propagators of the quasiparticles. Wavy lines stand for 
the charge bosons, and solid lines are the spin bosons. Con- 
densations of bosons are represented by lines terminated at 
crosses. Coupling constant g\ a g\ a ^ij is represented by open 
polygons. Here, we do not distinguish holons and doublons. 
Spins are also not distinguished in the diagram, for simplicity. 



(T-3) Diagrams that do not contain time dependence of 
quasiparticles but do violate the local conservation 
(FigEf). 

(T-4) Diagrams that neither include time dependence 
of quasiparticlcs nor violate the local conservation 
(FigEg). 

Here we present our guiding principle to take account 
of boson fluctuations: We exclude (T-l) because it vio- 
lates the local conservation when the quasiparticles dy- 
namically fluctuate. On the other hand we retain dia- 
grams belonging to the categories (T-2), (T-3), and (T- 
4). The reason to retain these diagram is as follows. 
The diagrams in the category (T-2) do not violate the 
local conservation, when bosons fluctuate. Therefore, we 
take the diagrams in this category into account. On the 



other hand, the diagrams in the category (T-3) do vio- 
late the local conservation. However, in these diagrams, 
quasiparticles enter as time averaged Green's functions. 
Therefore, quasiparticles feel the time averaged bosonic 
motions. The real violation of the local conservation oc- 
curs only when a dynamical quasiparticle process is in- 
duced by fluctuating boson hoppings. On the contrary, 
the real violation does not occur when the quasiparticles 
emerge as the time averaged quantities as in the case of 
(T-3). This is the reason to retain the diagrams in the 
category (T-3). Since (T-4) does not violate local con- 
servation, we retain it. 

For the slave-particle formalism of correlated fermion 
systems, it is well-known that fluctuations of gauge fields 
play an important role on reinforcing the local constraint 
imposed on slave particles^. It was pointed out by Joli- 
coeur and Lc Guillou that the Kotliar-Ruckenstein for- 
malism has the U(l)xU(l)xU(l) gauge symmetry^. It 
comes from the phase symmetry of the slave bosonic par- 
ticles, namely, e i , p ia , and d i . 

In our theory, we will treat fluctuations of such 
phases together with fluctuations of the amplitude of 
the condensation fraction of these slave particles, by 
using the Bogoliubov prescription. Therefore, the 
phase fluctuations are taken into account, although the 
U(l)xU(l)xU(l) gauge structure is not strictly con- 
served. 



(a) 



(b) 



FIG. 6: Diagrams for the transformed Lagrangians. Solid and 
dashed lines with arrows represent propagators of the quasi- 
particles and cofermions, respectively. Wavy lines stand for 
the charge bosons, and solid lines are the spin bosons, (a) The 
diagram represents the Lagrangian C' a i (Eq. (|A.I6|) h (b) The 
diagrams stand for terms in the Lagrangian C ' a 2 fEq. (|A.17|l V 



bosons depicted in FigJS^, 

£ a = ^C,L(r)(T ff ). y C ja (r) 



(A.I3) 



y<7 



as 

exp 
where 



jl[d±ld± ilT e-Io d ^ 



det [T T T ;J 



a — £ al + £ o2j 



-iA.u) 



(A.I5) 



b. Stratonovich- Hubbard tr as formation 



We introduce Grassmannian valuables (or fermionic 
fields) T jCT = (ip ia ,X ia ) T that stands for the cofermions 
as are discussed in the main article, by using a following 
identity 



Y[d±ld± ia e A = det 



(A.9) 



where matrices T„ are defined as 



PiaPja PiaP)a 
PzaPja PiaP)a 



(A.IO) 



and 



C'ai = £*J>)C r «0« t i*( r )' (A.16) 



C'a2 = -^{^(rKT^yf^fr) 



tl(r)T t3 C ]a (r)} 



(A.17) 



These transformed Lagrangian C a i (Fig|Bk) and C ' a 2 
(FigJ&b) lead to the cofermions' self-energy and the hy- 
bridization between the quasiparticles and cofermions, 
respectively after integrating out the fluctuating bosonic 
degrees of freedom. It will be discussed below, by using 
a set of the Dyson equations. 



ij<r 

x {T a )a(r ja {T)-c jtr (T) 



(A.11) 



Here we use vector notations as are defined in the main 
article as 



C 



c. Prescription for self- consistent procedure and Green's 
functions 

Here we construct approximated Green's functions for 
the Gaussian fluctuations of the bosons, quasiparticles, 
and cofermions by using a set of the Dyson equations as 
is depicted in FigJT] Thick lines and thick wavy lines 
stand for the Green's functions of the charge bosons 

L(eld t ) T . (A.12) A ab (r,r) = - (T/3?(r)/?f(0)) , and the spin bosons 



The identity Eq. ([A~9l) gives the transformation for 
a coupling term of the quasiparticles and fluctuating 



C ab (r,T) = - /T0?(r)<^ t (O)') , respectively, where a,b 



1,2, r = i - j, (PlPf) = (e^d,), and 
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{ViaiPto )■ Thick lines with arrows represent the quasi- 
particles (r, r) . On the other hand, Thin lines and 
thin wavy lines represent bare propagators of the charge 
bosons AQ b (r, r), the spin bosons CQ b (r, r), respectively, 
determined by £g , in which self-energy effects arc not 
taken into account. Thin lines with arrows stand for bare 
propagators of the quasiparticles Q^J (r, r) determined by 

A) = £ /t(r)[A^ + (oaUAfiAr), (A.18) 



where Coo- = gl a gl a (p 0cr eo + d p ow ) 2 . The Lagrangian C Q 
is obtained by decoupling the fluctuating bosons from 
the Lagrangian Co (Eq. (|A.6[) b Thick and thin dashed 
lines stand for the cofermions' propagators F ab and bare 
propagators Jq 6 = 5 a ^/e (e — > 0), respectively. 




FIG. 7: (Same figure as Fig.l in the main article) Diagrams 
for the Dyson equations. Solid and dashed lines with arrows 
represent propagators of the quasiparticles and cofermions, re- 
spectively. Thin wavy lines stand for the charge bosons, and 
thick solid lines are the spin bosons. Condensations of bosons 
are represented by lines terminated at crosses. Coupling con- 
stant giadia'tii i s represented by open polygons. Here, we 
do not distinguish holons and doublons. Spins are also not 
distinguished in the diagram, for simplicity. 



In the set of Dyson equations (Fig[7]), we neglect the 
coupling between charge and spin bosons described by 

propagators such as ^pt^Y at the Gaussian level, since 

these coupling terms are higher order contributions. Be- 
low we explain that the coupling gives higher order con- 
tribution with respect to hole-doping rate x, in proximity 
to Mott insulating states: Since operators including both 
charge and spin such as p\ (y e i do not conserve the electric 
charge, propagators such as (p\ cr e i ) should vanish in the 



Mott insulating phase, where the charge can not fluctu- 
ate. Therefore, the charge and spin excitations are well 
separated in the Mott insulating phase. 

When hole carriers are doped, p^e^ can have a non 
zero expectation value, at most, scaled by condensate 
fraction of holons eo which gives a rough estimate of the 
amplitude of charge fluctuations. From a relation oc x 
held in the KR theory for the hole-doped case, we obtain 

Pia^i ^ ^ Futhermore, there is an additional con- 
straint for the coupling terms such as ^pjo-e-^: they do 
not appear alone in calculations of physical quantities. 
To conserve charge and spin on average, (pl r7 e i \ appears 

with \ e\p i( j) hi pair, for example. Therefore, the contri- 



bution of the coupling between charge and spin bosons to 
physical quantities is scaled by (y/x) 2 . It concludes that 
the coupling between the charge and spin bosons gives 
contributions as a higher order in terms of x in physical 
quantities. 

By solving the set of the Dyson equations, we obtain 
the propagators for the quasiparticles and cofermions. 
Here the bosonic degrees of freedom are taken into ac- 
count in self-consistent fashion, through the cofermion 
self-energy S^(r, r), and the amplitude of hy- 

bridization between the quasiparticles and cofermions, 
each of which we detail below. 

The Lagrangian for the cofermions is given by 



£ cf = -£ii(T) 



,(cf), 



(r,r) T ja (r), (A.I9) 



where T icr = (4 , ia->Xi C r) i s a vector notation for the 
cofermions, as is defined in the main article, and r = i—j. 
The cofermion self-energy X^ cf ■* (r, r) is a 2 x 2 symmetric 



matrix. 



S (cf) 



E' 



E 22 



(A.20) 



The hybridization between the quasiparticles and 
cofermions is described by 

4yb = £[* !»A^v(t) + /X(T)A£t 3 v(rX£.21) 

where A^ = (A^A^). As a result, the effective La- 
grangian for the quasiparticles and cofermions, £ ff, is 
given as 



C 



eS 



A) 



-cf 



Aiyb- 



(A.22) 



When the charge gap is relatively small, E^ 



ii 



->22 



and A 



(VO 



A^ hold approximately. When the charge 



gap collapses, £ 



li 



-22 



and A 



MO 



A& 5 hold ex- 



actly. In our results, we employ the approximate rela- 



tions £ 



Hi 2 and A, 



A 



('i-) 



A[f. Then 
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a cofermion mode (ip ka + Xfc cr )/v / 2 hybridize with quasi- 
particles through the amplitude A(fc), which is depicted 
in FigJT] as closed (blue) circles, where k is a momen- 
tum. The inverse of cofermion propagator (namely, the 
cofermion self-energy) for (ip ka + x ka )j\/2 is given as 



— [E CT (fc, ie n ) + K(k, ie n )] = lk ie n -a k + 0(^.23) 

where fermionic Matsubara frequency. 

Then, the Green's function for the quasiparticles is 
given as 



Q [ J\k,ie r 



LU + i5 - Coaf-k + M 



A(kY 



uj + i5) = G ( J'(k,uj + i6) 

i-l 



7fe(w + iS) - a k 



(A.24) 



where e k is the Fourier transformation of £y , and [i is the 
chemical potential. Here we note that the weights of the 
two quasiparticle bands split by the zero surface defined 
by uj = at /'Yk are not the same in our theory. 

In our calculations, we define the doping rate x by 
using the quasiparticle Green's function as 



T 



(A.25) 



where T stands for temperature and N s is the number of 
sites. 

The Green's function for the electrons, instead of the 
quasiparticles, is given as 



Tc la ( T )cUO) 



-(t4(t)v(o)) 



x (TfMfUo) 



= (ri(r)i J5 (0))^(r), (A.26) 

where the bosonic and fermionic degrees of freedom are 
decoupled, because the resultant action in our theory 
does not contain the hybridization between bosons and 
fermions. The quasiparticle Green's function is defined, 
as in the previous sections, as 



TL(r)fU0) 



The bosonic part in Eq. (|A.26[) is given by 



(A.27) 



Tzl(r)z j<r (0)) ^ 9L9L{m(r)-p ia (T)} 

x [pj»-W])> (A.28) 



where we use vector notation as b i = (ej,^), p\ a = 
{PiaiPia)- Because we adopt the boson dynamics in 



u) (meV) 

-200 -100 100 200 

2.0 




-0.4 -0.2 0.2 0.4 

U)/t 



FIG. 8: Quasiparticle density of states (DOS). Solid curves 
show our result for DOS as function of ui for x — 0.05. Dashed 
curve shows DOS for the non-interacting case with the same 
band parameter t' = 0.25t and the doping x = 0.05. 



which charge and spin bosons are decoupled, this bosonic 
part of the Green's function is rewritten as 



Tz\ a {T)z ia (Q) 



T-T. 



gt9L{T[bo-pU\Po,-b ] 



+9UL {T[bi(r)-pU^o,-bM 
+9l9h{T[b -p ia (T)]\pl{T)-bo] 
+9l9l (nb'dr) ■ PMMM ■ b,(r)])(A.29) 

where b = (e ,d ) and p 0a = (p 0o .,p -). If we retain 
only the first and second lines of the right hand side of 
Eq. ([A.29[) . the electron Green's function is reduced to 
that already obtained in Refl3ll The contribution of the 
fourth line of the right hand side of Eq. (|A.29p is small 
compared with these from other lines in Eq. (|A.29[) , and 
we ignore the fourth term. 



2. Supplementary result for quasiparticle density 
of states 



Here we show supplementary results for the in- 
dependence of the quasiparticle density of states (DOS). 
The w-dependence of DOS (solid curves in FigH} shows 
significant asymmetry around w = compared to the 
DOS for the non-interacting case (the dashed curve in 
Fig |8|). 

This asymmetry of the DOS naturally explain the 
asymmetric tunneling spectra with respect to the sign 
of the bias observed in the scanning tunneling mi- 
croscopy (STM) measurements of the hole-underdopcd 
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cuptrate o 25 ' 35 . Although in the STM measurements, the 
asymmetry is observed up to several hundreds meV, the 
present result offers a possible origin of the asymmetry 
observed especially up to 100 meV. Our result is in sharp 
contrast to the previous work by Anderson and Ong22, 
in which the quasiparticle weights of the states for the 
added electrons above fj, and the states for the removed 
electrons below fi arc different from each other. They 
claimed that the quasiparticle weights inevitably show 
a step-like singularity at the Fermi level in the proxim- 



ity to Mott insulators. In our theory, such singularities 
at the Fermi level are not needed for occurrence of the 
asymmetric DOS. On the other hand, the recent study by 
Nieminen et al, shows that layers other than the Cu02 
layers in cuprates^ play a considerable role on the tun- 
neling spectra and cause asymmetric spectra. Such ef- 
fects of the layers other than the Cu02 layers, which are 
ignored in our theory, will enhance the asymmetry of the 
tunneling spectra. 
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